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Abstract 

We investigate strong-coupling effects on normal state properties of an ultracold Fermi gas. 
Within the framework of T-matrix approximation in terms of pairing fluctuations, we calculate the 
single-particle density of states (DOS), as well as the spectral weight, over the entire BCS-BEC 
crossover region above the superfluid phase transition temperature T c . Starting from the weak- 
coupling BCS regime, we show that the so-called pseudogap develops in DOS above T c , which 
becomes remarkable in the crossover region. The pseudogap structure continuously changes into a 
fully gapped one in the strong-coupling BEC regime, where the gap energy is directly related to the 
binding energy of tightly bound molecules. We determine the pseudogap temperature T* where 
the dip structure in DOS vanishes. The value of T* is shown to be very different from another 
characteristic temperature T** where a BCS-type double peak structure disappears in the spectral 
weight. While one finds T* > T** in the BCS regime, T** becomes higher than T* in the crossover 
region and BEC regime. Including this, we determine the pseudogap region in the phase diagram 
of ultracold Fermi gases. Our results would be useful in the search for the pseudogap region in 
ultracold 6 Li and 40 K Fermi gases. 

PACS numbers: 03.75.Ss,05.30.Fk,67.85.-d 
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I. INTRODUCTION 



Ultracold atomic Fermi gases provide unique opportunities to investigate the crossover 
from the Bardeen-Cooper-Shrieffer (BCS) type superfluids to the Bose-Einstein condensa- 
tion (BEC) of tightly bound molecules |2j, y, IJ] in a unified manner One 
of the key ingredients to achieve this BCS-BEC crossover in Fermi gases is a Feshbach 
resonance [9|, which allows one to tune the pairing interaction from the weak-coupling BCS 



n 



limit to the strong coupling BEC limit [id, [l2j, ll3j, Q. Since the BCS-BEC crossover 



is a fundamental many-body problem, it has recently attracted much attention, not only 
in cold atom physics, but also in various research fields, such as condensed matter physics 
and high energy physics. In particular, this system is expected to be helpful for further 
understanding of high-T c cuprates, which has been one of the most challenging problems in 
condensed matter physics 15]. 

In the under-doped regime of high-T c cuprates, the so-called pseudogap phenomenon 
has been extensively studied flol. [if] ]. In this phenomenon, the single-particle density of 
states (DOS) in the normal state exhibits a dip structure around the Fermi energy. The 
temperature at which the pseudogap appears is referred to as the pseudogap temperature 
T*, which is higher than the superconducting phase transition temperature T c . In the 
region between T* and T c , various anomalies have been observed in physical quantities, such 



as nuclear spin-lattice re 



axation rate (NMR-T^ x ) 17|, and angle-resolved photoemission 



spectroscopy (ARPES)[18|]. As the origin of the pseudogap, possib ility of pre 
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ormed pairs 



2J, UM- However, 



due to strong pairing fluctuations has been proposed 
because of the complexity of high-T c cuprates, other scenarios have been also discussed, 
such as antiferromagnetic spin fluctuations 



25 



26| and a hidden order 



27| . Thus, a simple 



system only having strong pairing fluctuations would be helpful to confirm whether or not 
preformed pairs are responsible for the pseudogap formation in high-T c cuprates. 

In this regard, the cold Fermi gas system meets this demand. This system is much cleaner 
and simpler than high-T c cuprates, and the pairing mechanism associated with a Feshbach 
resonance has been well understood. The BCS-BEC crossover is dominated by strong pairing 
fluctuations, so that one can focus on how they affect physical quantities. Indeed, effects 
of pairing fluctuations on single-particle _spectral weight have been theoretically studied by 



many researchers 



21 



22 



23|, 



24 



2a 



30, 



31 



321 ] . They clarified that pairing fluctuations 
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lead to a BCS-type double peak structure in the spectral weight above T c , which is a signature 
of pseudogap phenomenon. They also found that the two peaks in the spectral weight merge 
into a single peak at high temperatures. In Ref. 24J], detailed analysis on the spectral weight 
above T c has been carried out over the entire BCS-BEC crossover, and, in the BEC regime, 
the deviation from the BCS-type behaviors due to an asymmetric double peak structure 
has been pointed out. Since a photoemission-type experiment has recently become possible 



in cold atom physics 



331 ]. we can now examine strong-coupling effects on single-particle 



excitations within the current experimental technology. Although cold Fermi gases are 
not exactly the same as high-T c cuprates (e.g., pairing symmetry), the study of pseudogap 
phenomenon in cold Fermi gases is expected to be useful for further understanding of the 
underdoped regime of high-T c cuprates. 

In this paper, we investigate pseudogap behaviors of an ultracold Fermi gas above Tn. 
Including pairing fluctuations within the T-matrix approximation developed in Refs. 22|, |24| , 
we systematically examine how the pseudogap develops in DOS, as well as the spectral 
weight, over the entire BCS-BEC crossover region. We determine the pseudogap temperature 
T* at which the dip structure in DOS vanishes. We show that T* is quite different from the 
temperature T** where the double peak structure in the spectral weight disappears. In the 
BCS regime, we find that T* > T**. However, T** becomes higher than T* in the crossover 
region and BEC regime. Including this, we determine the pseudogap region in the BCS-BEC 
crossover phase diagram in terms of temperature and the strength of pairing interaction. 

This paper is organized as follows. In Sec. [Ill we explain our model and formulation 
to study pseudogap in DOS and spectral weight. In Sec. IHIl we examine the pseudogap 
structure in DOS. Here, we show how the pseudogapped DOS continuously changes into 
fully gapped one, as one passes through the BCS-BEC crossover region. We determine the 
pseudogap temperature T* from the temperature dependence of DOS. In Sec. HVl we examine 
strong-coupling effects on the spectral weight. We introduce another pseudogap temperature 
T** from the temperature dependence of spectral weight. We also discuss difference between 
T* and T**. Throughout this paper, we take h = k-Q = 1. 
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II. MODEL AND FORMALISM 



We consider a three-dimensional uniform Fermi gas, consisting of two atomic hyperfine 
states described by pseudospin o =f, |. So far, all the experiments on cold Fermigases are 
using a broad Feshbach resonance to tune the strength of a pairing interaction 
In this case, detailed Feshbach-induced pairing mechanism is known to be not crucial as 
far as we consider the interesting BCS-BEC crossover regime, and one can safely use the 
ordinary single-channel BCS model, described by the Hamiltonian, 

H =J2 ZP C U C P° ~ U C p+ q /2T C - P +q/2.l C -P'+<7/2!Cp'+<j/2t- (1) 

p,a q p t p' 

Here, c pcr is the annihilation operator of a Fermi atom with the pseudospin a and the kinetic 
energy £ p = e p — fx = p 2 /2m — //, measured from the chemical potential /i (where m is 
an atomic mass). —U (< 0) is an assumed tunable pairing interaction associated with a 
Feshbach resonance. It is related to the s-wave scattering length a s as 

4ira s U 



where u c is a high-energy cutoff. Since the strength of an interaction is usually measured 
in term of the scattering length a s in cold atom physics, Eq. (j2J) is useful in comparing 
theoretical results with experiments. In this scale, the weak-coupling BCS limit and strong- 
coupling BEC limit are characterized as (kpas)^ 1 <C — 1 and (kpas)^ 1 3> +1, respectively 
(where kp is the Fermi momentum). The region —1 < (/cpa s ) _1 < + 1 is referred to as the 
crossover region. The center of the crossover region ((/cf^s) -1 = 0) is called the unitarity 
limit 35]. 

To discuss strong-coupling effects in the BCS-BEC crossover regime above T c , we include 



22 



24|. Namely, we consider the 



pairing fluctuations within the T-matrix approximation 
single-particle thermal Green's function, 

G p (iu n ) = 1 , ; r, (3) 

where u n is the fermion Matsubara frequency. The self-energy correction E(p, iu n ) describes 
effects of pairing fluctuations, which is diagrammatically given by Fig. [T](a). In Fig. (TJ the 
solid lines are the free fermion propagator, 

G° p (iu n ) = (4) 



(a) 
(b) 
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FIG. 1: (a) Self-energy correction S(p, iuJ n ), and (b) particle-particle scattering matrix T(q, iv n ), in 
the T-matrix approximation. The solid and wavy lines represent the non-interacting Fermi Green's 
function Gp(iu> n ) and pairing interaction —U, respectively. 

Although this T-matrix theory does not treat the single-particle Green's function self- 
consistently, Ref. 24j has shown that it can correctly describe the smooth crossover from 
the BCS regime to the BEC regime. We briefly note that the self-consistent T-matrix ap- 
proximation (where the full Green's function G is used in stead of G° in evaluating the 



self-energy) has 
crossover region 



3een recently employed to study the spectral weight and rf-spectrum in the 



32] 



Summing up the diagrams in Fig. Q^a), we obtain 



E(p,iuj n ) = Ty^T{q,iv n )G 



i(v n -u} n )6 



(5) 



where v n is the boson Matsubara frequency. The particle-particle scattering matrix T(q, w n ), 
which describes fluctuations in the Cooper channel, is diagrammatically given by Fig. [H(b) . 
The expression is given by 

-U 



T(q,iu n 



1 — UTl(q, iv n 

47ra, 



ni 



rn 



U(q,iu n ) - 5^^- 



(6) 



In the last expression, the ultraviolet divergence coming from the contact pairing interaction 
has been absorbed into the scattering length a s 34|. U(q,iv n ) is the pair-propagator, given 
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by 



U(q,iu n ) = T^G p+q/2 (iv n + iu n )G _ p+q/2 (-iuj n ) 



p,u) n 



E 



1 ~ /(^p+q/2) ~ /(£p-q/2 



(7) 



where /(e) is the Fermi distribution function. 

To examine the pseudogap region, one needs to determine T c [jj fl, Q, 



for T r is obtained from the Thouless criterion 



241 ] . The equation 



36), T(g = 0,iz/ n = 0,T = T c 



1-1 



gives 



47ra s 



in 



E 



2(e p - /x) 



tanh — ^- 

2T 2 £p 



0, which 



(8) 



As pointed out by Nozieres and Schmitt-Rink|3], the chemical potential fi deviates from the 
Fermi energy ef in the BCS-BEC crossover. This strong-coupling effect can be conveniently 
included by solving Eq. (JSj) , together with the equation for the number N of Fermi atoms, 



N 



2T E 



GJiUr, 



(9) 



We show the self-consistent solutions of the coupled equations (JSj) and ([9]) in Fig. [2j 

In the normal phase above T c , we only solve the number equation to determine the 
temperature dependence of /i(T > T c ). The resulting fi(T) in Fig.[3]is used to calculate DOS 
p(oj), as well as the spectral weight A(p,u). They are obtained from the analytic continued 
Green's function as, respectively, 



— lm[G(p, iuj — > uj + iS)], 
n — ' 



(10) 



A(p, uj) = Im[G(p, iuj — > uj + iS)]. 



71 



(11) 



The analytic continued self-energy in G(p, iuj n — > uj + iS) has the form, 

E{p,u + i5)=X B + l-Yl I" dz nB ^ + f ^ q ^ 1m\T(q,iv n ^z + iS)], (12) 



z — {uj + iS) — £, 



q-p 



where Ub{£) is the Bose distribution function. Eh = — (U/2) J2 P fi^p) * s the Hartree term, 
and the last term in Eq. (JT2l) describes fluctuation correction to single-particle excitations. 

Before ending this section, we comment on the T-matrix theory used in this paper. In 
the BCS-BEC crossover literature, the so-called Gaussian fluctuation theory developed by 
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FIG. 2: Self-consistent solutions of the coupled equations (jSJ) and Q in the BCS-BEC crossover 
('TMA' in the figure), (a) phase transition temperature T c . (b) chemical potential n(T = T c ). In 
panel (b), \x is negative when (/cpa*,) -1 > 0.35. 'BCS' and 'NSR' are the weak-coupling BCS result 
and the NSR result, respectively. 

Nozieres and Schmitt-Rink (NSR)j^, 4] has been also used. The present T- matrix theory 
is a natural extension of this to include higher order pairing fluctuations. Indeed, the In- 
equation (jSJ) is common to the two theories, and the NSR number equation is also obtained 
from Eq. (jUJ), by expanding G p (iuj n ) in Eq. (Q up to 0(E), as 



Of K (iu n ) = G Jiou n ) + G°Muj n )J:(p,tuj n )G Muj n ). 



(13) 



The two theories essentially give the same BCS-BEC crossover behaviors of T c and /i(T = 
T c ), as shown in Fig. [2J In particular, both theories correctly describe the strong-coupling 
BEC limit, where the superfluid phase transition is dominated by BEC of N/2 tightly bound 
molecules (which leads to T c = 0.218T F [3]) and 2\fi\ equals the binding energy of a two-body 
bound state -Ebind = l/ma 2 s |2j. However, when one uses Gp SR (iu n — > to + i5) in calculating 
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FIG. 3: Calculated chemical potential /j, above T c in the BCS side (a) and BEC side (b). Each line 
starts from T c . We will use these results in calculating the density of states and spectral weight in 
Sees. Ill and IV. 

Eq. ([TO]) , unphysical results are obtained. The NSR theory overestimates the suppression 



of DOS around to — 0, leading to a negative DOS around u = in the crossover region 



37] 



The NSR theory also gives an unphysical divergence of DOS at ui = (although we do not 



explicitly show this in this paper) 37]. Thus, although the NSR theory can describe the 
BCS-BEC crossover behaviors of T c and /i, one needs to be careful in considering single- 
particle properties in the BCS-BEC crossover. Since this problem is absent in the present 
T-matrix theory, we employ this framework to examine DOS and the spectral weight in this 
paper. 
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FIG. 4: Density of states at T c . (a) BCS side ((fcpa,)" 1 < 0). (b) BEC side ((A^a*) -1 > 0). 

III. PSEUDOGAP IN SINGLE-PARTICLE DENSITY OF STATES 

In this section, we discuss the pseudogap phenomenon in DOS. Figure H] shows DOS in 
the BCS-BEC crossover at T c . Starting from the weak-coupling BCS regime, a pseudogap 
develops around uj = 0, as one increases the strength of the pairing interaction. Since the 
superfluid order parameter vanishes at T c , this dip structure purely originates from pairing 
fluctuations. 

The reason why the fluctuation correction described by the self-energy in Eq. Q causes 
the pseudogap in DOS can be easily understood by noting similarity between Eq. ([3]) and 



the Green's function in the mean-field BCS theory 



38] 



G P (^)--^ + a + A2 > ( 14 ) 

where A is the superfluid order parameter. Assuming that pairing fluctuations are strong 
around q = v n = (Note that T(q — 0, v n — 0) diverges at T c .), we may approximate 



Eq. © to 

E(p, iu n ) ~ S H - G° p (-^ n ) Aj g , (15) 

where Aj; g = — qVn [T(q, iis n ) + U]. Although G°_ p in Eq. (|T5|) does not involve the 
Hartree term Eh in the present T-matrix approximation, a better approximation would 
involve it in evaluating S. In this case, substituting Eq. (|T5|) into Eq. ([3]), we obtain 

where /i in £ p is replaced by /i + Eh- Since G Q _ p {—iuj p ) may be regarded as the hole Green's 
function, Eq. f|T6|) means that pairing fluctuations induce a particle-hole coupling. Compar- 
ing Eq. (fT6l) with Eq. ffT4l) . we find that A pg (which describes effects of pairing fluctuations) 
plays the same role as the BCS gap parameter A. Actually, dynamical effects of pairing 
fluctuations with q ^ and v n ^ smear the clear gap structure and coherence peak known 
in the mean- field BCS theory. However, in Fig. H](a), one can still see broad peaks around 
uj/e-f — ±0.2 (which correspond to the diverging coherence peaks at u = ±A in the BCS 
theory) when (/cpcis) -1 ^ — 0.4. Although the above discussion simplifies the treatment of 
pairing fluctuations, it would be helpful in understanding the reason why pairing fluctuations 
give the pseudogap structure above T c . 

While the pseudogapped DOS is very remarkable in the unitarity limit, it continuously 
changes into a fully gapped one in the strong-coupling BEC regime, as shown in Fig. H](b) . 
In the BEC regime where /i is negative ((kpag)^ 1 > 0.35), when we only retain the negative 
fi and ignore other strong-coupling effects, the DOS has a finite energy gap as 

i (<*>< M)> , 17 x 

= { ,/ 2 , ( 17 ) 



in 



In the BEC limit, 2|/i] equals the binding energy -Ebind = 1/™^ of a two-body bound state, 
which means that the energy gap in Eq. (fTTl) is directly related to the molecular dissociation 
energy. Since the intensity of DOS is almost absent below uj/e-f ~ 1.4 when (/cpcts) -1 = +0.8 
in Fig. BJb), the region of (/cf^s) -1 > 0.8 is considered to be close to an N/2 molecular gas, 
rather than an N atomic Fermi gas. 

However, we note that p(u < 0) still has small but finite intensity even when 
(fcpOs) -1 ^ 1.0, as shown in Fig. BJb), which means the existence of hole-type excitations. 
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The finite DOS in the negative energy region is absent when we ignore all fluctuation ef- 
fects except for the negative \x (See Eq. (fTTI) .). Since the concept of hole is characteristic 
of many-fermion system, one finds that, although the BEC region around (/cp^s) 1 — +1 
is dominated by two-body molecular bosons, the character of many-fermion system still re- 
mains to some extent there, leading to the finite p(u < 0). We also find this by simply 
employing Eq. (Ti~6[) to calculate DOS in the BEC regime (// < 0), which gives 

yj\l» 2 -Xk-\n\ (" < ~ V a p S + H a ) 



p{u) = { 




^ 2V2n 2 y/u' 2 - A 





A 2 

pg 


UJ 




A 2 

pg 





When the two-body binding energy -Ebind — l/nia 2 (~ 2|/x|) is much larger than the 'charac- 
teristic energy' A pg , one may ignore A pg in Eq. ([TBI) . In this extreme BEC limit, the upper 
branch in Eq. ([18]) reduces to Eq. ffTTj) . and the lower one vanishes, as expected. 

Figures [5] and [6] show DOS above T c . The pseudogap structure in DOS becomes obscure 
at high temperatures due to weak pairing fluctuations. The dip structure eventually vanishes 



at a certain temperature, which we define as the pseudogap temperature T* [40 1 . 

Figure [7J shows the resulting pseudogap temperature T* in the BCS-BEC crossover. Start- 
ing from the weak-coupling BCS regime, T* monotonically increases. However, T* is still 
lower than T c calculated in the mean-field BCS theory ('BCS' in Fig. [7J). Although the mean- 
field T c is sometimes considered as a characteristic temperature where preformed pairs are 
formed, our result shows that the pseudogap actually starts to develop in DOS from lower 
temperature. 

We note that, although the fact that the pseudogap disappears at T* is common to the 
entire BCS-BEC crossover region, the detailed way of disappearance is somehow different 
in between the BCS regime and crossover-BEC regime. In Fig. [5(a), the pseudogap around 
uj = is simply filled up at high temperatures. The shape of DOS then becomes close to 
DOS of a free Fermi gas, 

pH = ^^V^ + P (w > -//). (19) 

Namely, as far as we consider DOS, the system may be regarded as a (weakly interacting) 
normal Fermi gas above T*. On the other hand, in the BEC side shown in Fig. [61 in 
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FIG. 5: Temperature dependence of the density of states p(co) in the BCS side. T c in each panel 
equals (a) 0.112ep, (b) 0.146£f, (c) 0.183£f, and (d) 0.217£f- In this figure and Figj6l we have 
offset the results for T > T c . The short horizontal line near each result is at p{oj) = 0. 

addition to the enhancement of DOS around uj = 0, the lower peak is suppressed at high 
temperatures. In the unitarity limit (Fig. EJ^a)), when the pseudogap is completely filled up, 
DOS still has a different shape from DOS of a free Fermi gas. In the BEC regime where 
fi < 0, Figs. 0(c) and (d) show that DOS above T* has a finite intensity in the negative 
energy region, in contrast to Eq. ffTT|) . These results indicate that pairing fluctuations still 
affect single-particle excitations above T* in the BEC side, although the depression of DOS 
around u = is absent. Indeed, in Sec. IIVt we will show an evidence of such fluctuation 
effects in the spectral weight in this regime. 
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FIG. 6: Temperature dependence of DOS in the BEC side. T c in each panel equals (a) 0.244£p, 
(b) 0.259e F , (c) 0.262e F , and (d) 0.255e F . 



IV. PSEUDOGAP IN SPECTRAL WEIGHT 



It has been pointed out 



21 



22 



23 



24j that pairing fluctuations cause a BCS-type double 



peak structure in the single-particle spectral weight A(p,uj). In this section, we examine 
how this strong-coupling effect is related to the pseudogap in DOS discussed in Sec. IIHI 

Figure [8] shows the intensity of spectral weight A(p,uj) at T c in the energy-momentum 
plane. In the BCS side (panel (a)), in addition to the particle branch at uj ~ £ p , we can see a 
weak peak line of a hole branch at uj ~ — £ p . The intensity of the particle branch is suppressed 
around uj = 0, where it intersects with the hole branch and the level repulsion between them 
occurs. The resulting structure is similar to the BCS spectral weight [2l|, |22|, |23|, |24|, |41| . 
given by[38j 

^bcs(p, oj) = u 2 p 5(uj - E p ) + v 2 p 6(uj + E p ), (20) 



where u% = (1 + f p /£ p )/2, v 2 p = (1 - £ p /E p )/2, and E p = ^ + A 2 is the Bogoliub 
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FIG. 7: Pseudogap temperature T* determined from DOS in the BCS-BEC crossover. We also 

weight 



391]. T* 



plot another pseudogap temperature T** where the double peak structure in the spectra 
vanishes. 'BCS' is T c = (87/vre 2 )e F e 7r/2fcFas in the mean-field BCS theory (where 7 = 1.78) 
or T** gives the boundary between the pseudogap regime (PG) and normal Fermi gas regime (NF). 
2|/i| (~ -Ebind) hi the BEC regime gives the characteristic temperature below which thermal dis- 
sociation of bound molecules are suppressed. Namely, T ~ 2|//| physically describes the boundary 
between PG and molecular Bose gas regime (MB). 

quasiparticle excitation spectrum. For a given momentum p, Abcs(p, has two peaks at 
uj = ±E P . The negative energy branch at 00 — —E p given by the second term in Eq. fl20|) 
is dominant in the low momentum region p <C p-p (where u p <C v p ). On the other hand, 
the positive energy branch (u = +E P ) becomes crucial when p 3> kp (where u p ^> v p ). 
The existence of two branches can be understood from the Bogoliubov transformation c p | = 
u plp1 + v p1^- p i {ipc is an annihilation operator of a quasiparticle with momentum p and 
spin cr), which indicates that the annihilation of an atom is accompanied by creation and 
annihilation of Bogoliubov excitations [42]. The minimum energy gap 2 A between the two 
branches uj = ±E p is obtained at the Fermi level p = k-p. Since the simplified Green's 
function in Eq. (fT6ll has the same form as Eq. (TT41) . Eq. (fT6l) gives rise to the spectral weight 
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FIG. 8: Calculated intensity of the spectral weight A(p, oS) at T c in the energy-momentum plane, (a) 



BCS side ((kpa. 



\-l 



-0.6). (b) Unitarity limit ((/cf^)" 1 = 0.01). (c) BEC side {(k F a s 



0.6). 



similar to the BCS type in Eq. (12"0|) . where the superfluid gap A is now replaced by the 
pseudogap A pg , describing effects of pairing fluctuation. The minimum value 2A pg of the 
pseudogap energy is obtained at p ~ hp in this case. From this reason, the double peak 
structure in Fig . [8](a) is found to come from the particle- hole coupling due to strong pairing 
fluctuations 24]. In addition, they also induce finite lifetime of quasiparticle excitations, 
leading to finite widths of the two peaks in A(p,u) 24|. This feature is absent in the BCS 
spectral weight in Eq. fl2"0~l) . which has two 5-functional peaks at u — ±E P . As a result, 
A(p,u) at the momentum where the minimum peak-to-peak energy is obtained has finite 
spectral weight between the two peaks, as shown in Fig. [9j giving finite intensity of DOS 
inside the pseudogap. This gapless double peak structure is referred to as the pseudogap in 



the spectral weight in the literature 21 



22 



23 



2i|. 



This pseudogap structure in the spectral weight becomes remarkable, as one approaches 
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FIG. 9: Spectral weight A(p, to) as a function of uj. In each panel, we take the momentum where 
the peak-to-peak energy becomes minimum: (a) p/kp = 0.91, (b) 0.83, and (c) 0.01. 

the unitarity limit. In this limit, strong pairing fluctuations also broaden the spectral peaks, 
as shown in Fig. E^b). In the BEC regime (Fig. EJ^c)), the peak width of the upper branch 
shrinks. This is because the BEC regime is well described by a gas of tightly bound molecules, 
so that the upper branch simply describes their dissociation. Since the molecular formation 
simply occurs within two-body physics in the BEC limit, the peak of the lower branch (which 
is an evidence of many-body physics) is low and broad in Fig. E^c). 

These different behaviors of upper and lower peaks in the BEC regime can be directly 
understood from the imaginary part of the self-energy correction. Using the fact that the 
particle-particle scattering matrix T reduces to the Bose Green's function in the BEC limit 



as 



24J 



S7T 1 . 

r( 9 ,M/ n = — -g (21) 

m z a s iv n - E° 

(where = q 2 /Am — \xb is the energy of a molecule measured from the molecular chemical 
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potential \xb — 2/i + l/(ma 2 ) ~ 0), we can approximately evaluate the imaginary part of 
the analytic continued self-energy in Eq. (|12p as 

■2 



lmZ(p,oj + i5) = —J2n B (E*)5 (u - (E* - £ q „ p )) 

Q 



AT , 
In 



a s p 



2 

x ^+wth-w), (22) 

where Au; = u; th — cj, and u th = fi — fi B ~ — \/2ma 2 . Since Im£(p, uj + i5) directly gives 
the peak width of the spectral weight, the first line in Eq. (1221) indicates that, in the BEC 
regime, the peak widths are dominated by molecules excited thermally with finite center of 
mass momentum q ^ 0. Since Eq. fl22|) vanishes when oj > p 2 /2m + uJth — p 2 /2m — l/2ma 2 , 
the upper branch around uj = £ p (> 0) appears as a sharp delta-function peak in the spectral 
weight in the BEC limit. This is consistent with the sharp upper peak in Fig. [9]^c). 

On the other hand, expanding Eq. ( 1221) around the lower branch, u = £ p , one obtains 

4T , / m „ 2 



Im£(p, uj + iS) ~ — In -r^Su 2 , (23) 

where 8uj = uj — (— £ p ). Equation (|23|) shows that the imaginary part of the self-energy 
logarithmically diverges along the lower branch uj = — £ p . Thus, the lower peak is smeared 
out in the BEC limit. 

As one increases the temperature, Fig. [9] shows that the double peak structure gradually 
becomes obscure to eventually vanish at a certain temperature (= T**). Regarding T** as 
another pseudogap temperature [43j, one might expect that it is deeply related to T* defined 
from DOS, because DOS is given by the momentum summation of the spectral weight. 
However, when we compare T** with T* in the BCS-BEC crossover, they are very different 



from each other, as shown in Fig. [71 While one sees T* > T** in the BCS side2j|, T** 
becomes higher than T* in the BEC side ((/cp a s) _1 ^ _ 0.07). 

In the BCS side, when T > T**, since pairing fluctuations are still strong near the Fermi 



surface, the single peak in the spectral weight at p ~ y / 2m/I is broad and the peak height is 
low, compared with the cases of higher and lower momenta, as shown in Fig. [TUJ This low 



peak height at p ~ v2m/i directly affects the density of states around uj = 0, leading to the 
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FIG. 10: (a) Intensity of the spectral weight A(p,u>) in the BCS side ((^f^s) -1 = —0.6). We set 
T/T c = 1.03, at which the dip structure can be clearly seen in DOS. (b) A(p, uS) as a function of 
uj. The momentum p is taken to be p/k-p = 0.91 (solid line), 0.83 (dashed line), and 0.97 (dotted 
line). 

dip or pseudogap structure in p(u) in the region T** < T < T*. We briefly note that the 
result of T* > T** in the BCS side agrees with the previous work [2 11] . 

On the other hand, although the double peak structure still exists when T > T* in the 
BEC side, the intensity of the lower peak is very weak and broad (See Fig. [TTJ), because the 
system is close to a gas of two-body bound molecules. Thus, the existence of lower peak is 
easily smeared out in the momentum summation in calculating DOS, p(oj) = ^2 p A(p,oj). 

To see the physical backgrounds of T* and T**, it is convenient to recall that, when pairs 
are formed above T c , the lifetime of Fermi excitations becomes short due to strong tendency 
to form pairs, leading to a broad quasi-particle peak in the spectral weight A(p,u). In 
addition, preformed pairs also induce the particle-hole coupling, which gives the double 
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FIG. 11: (a) Intensity of spectral weight A(p,uj) in the BEC side ((/cfo s ) 1 = +0.4). We take 
T/T c = 1.53, at which the pseudogap structure is absent in DOS. (b) A(p, lu) as a function of u at 
p = O.Olkp. 

peak structure in A(p, lu). Between the two effects associated with pair formation, while T** 
is directly related to the latter by definition, the former is crucial for T*: In the BCS regime, 
since the peak-to-peak energy in A(p,cu) is small, the double-peak pseudogap structure is 
easily smeared out by the lifetime effect, namely, the broadening of two peaks. On the 
other hand, DOS around lu = is suppressed, when the height of quasiparticle peak at 
lu ~ is lowered by the broadening effect. As a result, one obtains T* > T** in the BCS 
regime, and one may use T* as the characteristic temperature where preformed pairs are 
formed. The double peak structure can be clearly seen in A(p,u) in the crossover-BEC 
regime, because the peak-to-peak energy becomes larger than the peak widths. However, as 
discussed previously, the lower peak becomes very broad and the weight becomes small in the 
BEC regime, reflecting that the system is close to a gas of two-body bound molecules. Thus, 
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one cannot see the dip structure in DOS even when the particle-hole coupling induce the 
double peak structure in A(p,u) below T**. As one further decreases the temperature, the 
lower peak in A(p, u) shrinks and the peak height increases, because the system approaches 
the superfluid phase. This clearly enhances the intensity of DOS in the negative energy 
region, leading to the dip structure below T*(< T**). 

The different behaviors of two pseudogap temperatures T* and T** imply that the pseu- 
dogap region may depend on what we measure. When we consider a quantity where DOS 
is crucial, T* would give the boundary between the pseudogap region and normal Fermi 
gas regime. On the other hand, when we consider a quantity dominated by the spectral 
weight, T** would be observed as the boundary between the two regions. While the specific 



33[ is 



heat is an example of the former quantity, the recent photoemission-type experiment 
considered to be a latter example. 

We note that, when the temperature is lower than the binding energy E^md — 
of a two-body bound molecule in the BEC regime, thermal dissociation of molecules is 
suppressed. In this sense, one may regard this regime as a molecular Bose gas, rather than 
a (strongly-correlated) Fermi gas. Including this, we obtain the phase diagram in Fig. UJ In 
this figure, the pseudogap regime is the region surrounded by T* or T**, T c and 2|//|. We 
briefly note that except for T c , other temperatures T*, T**, and T = 2|/x|, are all crossover 
temperatures without accompanied by any phase transition. 



V. SUMMARY 



To summarize, we have investigated the pseudogap behaviors of an ultracold Fermi gas 
in the BCS-BEC crossover above T c . We have calculated the single-particle density of states 
(DOS), as well as the single-particle spectral weight, including pair fluctuations within 
the framework of T-matrix approximation. We showed how the pseudogap structure ap- 
pears/disappears in DOS above T c in the BCS-BEC crossover region. Starting from the 
weak-coupling BCS regime, while the pseudogap in DOS becomes remarkable near the uni- 
tarity limit, it continuously changes into a fully gapped DOS in the BEC regime. 

We determined the pseudogap temperature T* as the temperature when the dip structure 
in DOS disappears. We also introduced another pseudogap temperature T** at which the 
double peak structure in the spectral weight vanishes. We showed that, although both the dip 
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structures in DOS and the double peak structure in the spectral weight originate from pairing 
fluctuations, their values are very different from each other in the BCS-BEC crossover. 
While one finds T* > T** in the BCS side ((/c F a s ) _1 ^ 0), T** becomes much higher than 
T* in the BEC side ((/cp^s) -1 ^ 0). This means that the pseudogap region may depend 
on the physical quan tities which we measure. In particular, since the recent photoemission- 
type experiment [331] is related to the spectral weight, one expects that T** would work 
as the pseudogap temperature in this experiment. Including T* and T**, we determined 
the pseudogap region in the BCS-BEC phase diagram with respect to temperature and 
the strength of pairing interaction. Since the pseudogap effects are crucial in understanding 
strong-coupling Fermi superfluids, our results would be useful in the search for the pseudogap 
region in the BCS-BEC crossover regime of ultracold Fermi gases. 
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